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ABSTRACT 

The conditions under which high frequency (f»f u h) long-range extraordinary-mode discrete 
field-aligned echoes observed by the Radio Plasma Imager (RPI) on board the Imager for 
Magnetopause-to-Aurora Global Exploration (IMAGE) satellite in the plasmasphere are 
investigated by ray tracing modeling. Field-aligned discrete echoes are most commonly observed 
by RPI in the plasmasphere although they are also observed over the polar cap region. The 
plasmasphere field-aligned echoes appearing as multiple echo traces at different virtual ranges 
are attributed to signals reflected successively between conjugate hemispheres that propagate 
along or nearly along closed geomagnetic field lines. The ray tracing simulations show that field- 
aligned ducts with as little as 1% density perturbations (depletions) and less than 10 wavelengths 
wide can guide nearly field-aligned propagating high frequency X mode waves. Effective 
guidance of wave at a given frequency and wave normal angle OP) depends on the cross-field 
density scale of the duct, such that ducts with stronger density depletions need to be wider in 
order to maintain the same gradient of refractive index across the magnetic field. While signal 
guidance by field aligned density gradient without ducting is possible only over the polar region, 
conjugate field-aligned echoes that have traversed through the equatorial region are most likely 
guided by ducting. 

1. INTRODUCTION 

It has long been shown that radio signals launched into ionospheric magnetic field-aligned 
density irregularities ( FAT) can be guided or ducted along magnetic field lines [e.g., Muldrew, 
1963; 1967; Platt and Dyson, 1989; James, 2000;], analogous to whistler mode waves being 
guided along magnetic field lines in the magnetosphere (e.g., see Smith [1961], Helliwell 
[1965]). Ducting of freely propagating electromagnetic waves has also been expected to occur in 
the magnetosphere [Calvert, 1981; 1995]. In the case of ducting, FAI can act like a waveguide 
that traps and channels wave energy along the background magnetic field. For high-frequency 
electromagnetic waves (f »f u h, where fuh is the upper hybrid resonance frequency), a region of 
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density depletion extending some distance along the magnetic field can trap wave energy by total 
internal reflections along a field-aligned density duct, as depicted in Figure 1. 

In analyzing data taken by the Radio Plasma Imager (RPI) on board the Imager for 
Magnetopause-to-Aurora Global Exploration (IMAGE) satellite, Reinisch et al., [2001] 
suggested another possible signal guiding mechanism due to the presence of field-aligned density 
gradients, such as might be present over the polar cap regions. In this scenario, field-aligned 
propagating radio signals can naturally follow the magnetic field line without the trapping action 
of an electron density duct. 

The RPI, the first instrument to perform long-range magnetospheric sounding successfully 
[ Reinisch et al., 2000], transmits coded electromagnetic waves over a selectable frequency range 
from 3 kHz to 3 MHz from a long dipole antenna. The far-field radiation pattern of each 
transmitted pulse covers nearly all directions. From a sweep over a series of sounding 
frequencies RPI often observes return signals with discrete virtual-range and frequency 
characteristics over widespread ranges of altitudes and magnetic latitudes varying from the polar 
cap to the equator [Fung et al., 2002; 2003]. Figure 2a shows an RPI plasmagram (virtual range 
versus frequency) when the IMAGE spacecraft was just inside the plasmasphere. The virtual 
range is given by ct/2 where c is the free-space speed of light and t is the echo delay time. The 
echo intensity is color-coded. Three basic types of RPI echoes (field-aligned, diffuse, and 
resonance) are all recorded in Figure 2a. Plasma resonances observed by RPI and their 
usefulness in deriving local plasma and magnetic field strengths [e.g., Fung et al., 2003] have 
been discussed extensively by Benson et al. [2003], The diffuse echoes typically show range or 
frequency spreading whereas field-aligned echoes tend to appear in discrete traces in 
plasmagrams and in many cases multiple traces over the same frequency range. Based on 
direction analysis and the virtual range determination the ray paths of the diffuse and field- 
aligned echoes of Figure 2a are shown schematically in Figure 2b. Carpenter et al., [2002] and 
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Fung et al [2003] attributed the range spreading of the diffuse echoes to coherent (aspect- 
sensitive) scattering from field-aligned electron density irregularities ( FA1 ) having < 10% 
variations of the background density that range in size from 200 m to a few km as illustrated in 
Figure 2c, although Carpenter et al [2002] also pointed out the possibilities of refraction into and 
partial propagation along ducts having cross-field scales of several signal wavelengths. 

As Fung et al [2003] pointed out, the diffuse echoes seen from outside the plasmapause and the 
field-aligned guided echoes seen inside the plasmapause represent two manifestations of the 
presence of FAI. Although it is not clear what processes can produce echo-guiding conditions 
that can extend from hemisphere to hemisphere, the widespread occurrences of field-aligned 
guided echoes clearly point to the existence of such conditions in the diverse environments of the 
magnetosphere and ionosphere. 

Fung et al. [2003] showed that multiple traces of RPI discrete echoes of the type shown in Figure 
2a, particularly the so-called epsilon signature, are most likely due to ducted signals from 
conjugate hemispheres along closed field-aligned paths, analogous.to the mechanism generally 
invoked to explain similar ionospheric observations. Results from modeling of ducted signals in 
the ionosphere have been reported by Muldrew [1963; 1967; 1980] and by Platt and Dyson 
[1989] and analysis of direct bi-static observations have been performed by James [2000]. 

The purpose of this paper is to investigate the plasma conditions and wave characteristics 
necessary for high-frequency extraordinary mode waves (/»/*) to be guided along the 
geomagnetic field over large magnetospheric distances. In order to determine the circumstances 
responsible for guiding long-range discrete echoes observed by the RPI, we have performed ray- 
tracing calculations within a model magnetosphere. Results will be compared to studies of 
ducting in the ionosphere [e.g., James, 2000]. Aside from the similar dependence of the ducting 
mechanism on the cross-field wave refractive index profile, much of the results presented here 
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are not directly applicable to ducting of whistler mode waves, we will not attempt to compare 
our results with those of whistler ducts. 


2. RAY TRACING 

The computer code used for the ray tracing calculations in this study was originally developed in 
two dimensions by Shawhan [1967] and extended to three dimensions by Green [1984], The 
code simultaneously solves the set of first-order differential equations for ray propagation in 
spherical coordinates developed by Haselgrove [1955] and Haselgrove and Haselgrove [I960]. 
The ray tracing technique uses the cold plasma theory [Stix, 1992] to determine the index of 
refraction and compute the ray paths by using geometric optics formulations. Light rays in the 
freely-propagating L-O and R-X modes are reflected when they encounter their respective wave 
cutoffs, 

fLO = f P = 9jW e ( 1 ) 


Irx - 


L 

2 




( 2 ) 


where f p ,f g , and N e are the electron plasma frequency, gyro frequency, and electron density, 

respectively. Thus the frequency of the reflected signals can be used to determine the local 
electron density at the cutoff. 


To initiate a ray-tracing run, the code requires inputting a location of the source region (e.g., 
satellite position), wave frequency, the mode of propagation (. L-O or R-X) and the initial wave 
normal angle (angle between the wave vector A: and the local magnetic field B ) as initial 
conditions. Ray paths will then be calculated based on the propagation characteristics of a given 
mode in the model magnetosphere. Prior to the launch of the IMAGE spacecraft this code was 
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used extensively to predict the traces of free space mode echoes [Fung and Green, 1996; Green 
et al., 2000], According to Haselgrove [1955], the ray equations that describe the propagation of 
an electromagnetic wave in cold magnetized plasma are given by 
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where n 2 = Af t-t=n 2 r +nj+nj is the square of the index of refraction n with t=k r r+ k e 0+ k t tp 
being the wave normal vector. 


In order to simulate guided wave propagation in the plasmasphere, we have adopted a linear 
combination of several analytical models to prescribe the spatial variations of n. These include a 
diffusive equilibrium model (how the plasmas distributed) by Angerami and Thomas [1964], an 
ionosphere and plasmasphere compositional model (% of electrons, protons, He + , and O 4 ) by 
Kimura [1966], a plasmapause distribution (location and magnitude of the density gradient) by 
Aikyo and Ondoh [1971], a field-aligned density depleted duct (described below), and a dipole 
magnetic field model. These models provide the background plasma and structure of the 
plasmasphere over which our ray tracing simulations are performed. 
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In order to simulate the effects of a field-aligned density duct on the RPI transmitted signals, a 
duct density model is added. The total electron density inside a field-aligned density duct is 
modeled by 


N(r,0)=N DE \ 


1 + expi 


N, 


DE 


1 r-Lcos 2 O' 2 


(9) 


where Nde is the diffusive equilibrium model density, 5NNde is the fractional density change at 
duct center, r is radial distance of a field location, a is the “width” of the density duct, and 0is 
the magnetic latitude of the field location at r. It is clear from equation (6) that the density 
perturbation maximizes at r = L cos 2 8, i.e., along a field line identified by the dipole L value, 
corresponding to a maximum change in N e (x) across the background field (see Figure 1). For 
freely propagating modes, density depletions ( SN < 0) in the duct lead to enhancements of the 
refractive index «(x) at the duct center relative to the edges, such that total internal reflection 
within the duct is possible in accordance to Snell’s law [Born and Wolf, 1999]. On the other hand 
for whistler propagation, refractive index enhancements are caused by density increases {SN > 

0), although whistler ray directions are already confined generally to within 20° of the 
background magnetic field lines [Storey, 1953; Stix, 1992], 

3. RESULTS OF RAY TRACING CALCULATIONS 

As discussed above, RPI generates phase-encoded emissions at all wave normal angles at each 
sounding frequency. Using the models described above, ray paths of extraordinary-mode waves 
launched from a given position having different initial wave normal angles are calculated by 
solving equations (3) - (8). For simplicity, waves are launched from the field line at the center of 
the field-aligned density duct given by (9). Figure 3 shows 200-kHz X-mode rays (blue and 
green traces) launched in the magnetic meridian plane with initial wave normal angles 
ys — 0° - 360°, each separated from one another by 10°. The rays paths are launched from a 
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point at 20° magnetic latitude and L = 3.5. The duct model density has a depletion 8AW D e = - 
15% and a = 20A^or 30 km). 

Figure 3 shows that most rays with large wave normal angles escape the density duct and 
propagate into the plasmasphere. Rays that are launched obliquely to the magnetic field toward 
the interior of the plasmasphere are refracted and reflected at the right hand-extraordinary (R-X) 
cutoff surface, given by equation (2) and the red contours. Rays launched outward simply escape 
the plasmasphere. Only rays with relatively small wave normal angles (shown as green traces) 
are trapped and guided by the model density duct. Figure 4a shows a close-up view of nine 120- 
kHz rays launched at different wave normal angles from 0° to 72° at every 8° interval toward the 
northern hemisphere from the magnetic equator. Only rays with relatively small wave normal 
angles (y/ < 30°) are guided along the duct by total internal reflection and that only the ducted 
rays (green traces) can propagate to reach both conjugate hemispheres and are reflected back to 
the launching point, as shown in Figure 4b. 

3.1 EFFECT OF DUCT WIDTH 

In this section, we investigate how the effectiveness of ducting may depend on duct width, a. For 
a given field-aligned density duct with a fixed depletion SACVde < 0, varying the duct width 
means changing the rate at which density increases from the center toward the “edge” of the duct 
(\r-Lcos 2 » a), where the background density level JV DE is recovered (see equation 9). For a 
given frequency wave propagating with a small but finite wave normal angle, the width of the 
duct also governs how quickly the wave encounters significant density effects toward the edge of 
the duct. Wave trapping (and guiding) by a field-aligned duct, as depicted in Figure 1, occurs 
when a wave propagating obliquely near the center of the duct is refracted so that the wave 
propagates toward the edge of the duct with an angle of incidence <|> = |90 °-y\ greater than the 
critical angle <j) c for total internal reflection [Born and Wolf, 1999]. 
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If the wave frequency far exceeds the local cutoff frequencies (equations 1 and 2), however, the 
wave will simply go through the ducting region and escape. On the other hand, waves with 
frequencies below the cutoff frequencies are evanescent and thus cannot propagate. Therefore a 
field-aligned density duct of a given depletion and width can only trap waves in limited 
frequency range, 1 <f/fxx < P- Since the density depletion is typically of order of a few percent, 
we can normalize the wave frequency /by the local R-X cutoff frequency (equation 2) based on 
Nue and f g . Similarly we can express duct width as <r=TjA,or a = (ol X )=tj with X being the free- 
space wavelength, and determine p locally by launching rays at different positions along a duct. 

Figure 5 shows the maximum ducted frequency P =fmax/fyx of waves at different y / as a function 
of normalized duct width r| for two different depletion levels. The figure indicates that higher 
frequency waves are more easily ducted when they are closer to field-aligned propagation and in 
deeper depletion ducts. For a fixed depletion, a gentle transverse density gradient associated with 
a wider duct thus can only trap lower frequency waves because high frequency waves with 
shorter wavelengths can escape the duct more easily. For all duct sizes (with o > few X), all 
obliquely propagating waves trapped, say y/> 10°, can only have frequencies slightly greater 
than the local /oc- Figure 6 shows the ray-tracing results of P as a function of T| for different 
depletion levels when the initial wave normal angle is 0°, indicating the maximum ducted 
frequency ranges for different sets of duct parameters: 8MVde and t). 

As shown in Figure 5, waves at frequencies / < 10/tx and have different wave normal angles can 
be guided by field aligned density ducts with 8 AWde < 10%. The echoes having slightly different 
directions of propagation will have different time delays, causing a small spread in virtual range 
or a finite width of the discrete echo trace. In order to understand the virtual range spreading 
effect due to different signal propagation directions, we have performed the ray tracing 
calculations shown in Figure 7. Assuming a duct with r\ = 50 and 5 AW D e< 5%, Figure 7 shows 
that a 0.2 R E variation in virtual range can result from echoes (f/fnx =1.1) having different initial 
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wave normal angles (0° < yr< 1 8°)> enough to explain most typical widths of discrete (guided) 
echo traces observed by RPI in the magnetosphere [Fung et al., 2003]. Larger range spreading (~ 
0.5 Re), such as reported by Carpenter et al. [2002] and Fung et al [2003], may be caused by 
coherent scattering processes. 

3.2 SIMULATED PLASMAGRAM 

We can use ray-tracing modeling to construct simulated plasmagrams and show that observed 
plasmagram traces may in fact result from ducting. Figure 8 shows the results of ray tracing 
simulated plasmagrams when ducted rays are launched at the equator (panel A) and at 20° 
magnetic latitude (panel B). With only a few percent density depletion (8 AWde ~ 5%) for the 
entire density duct, waves at frequencies that exceed slightly the local R-X cutoff (equation 2) 
have relatively slow group speeds and thus have relatively long delay times or virtual ranges, 
causing the nearly-vertical portion of the trace as shown in panel A. Higher frequency waves 
with group speeds approaching light speed can reach deeper in the ionosphere where the density 
gradients are sharper, leading to the nearly flat or asymptotic trace appearance. 

From a source at 20° magnetic latitude within the model plasmasphere two distinct echo traces 
are observed in Figure 8B. The bottom trace of Figure 8B are simulated echoes from along the 
field line in the hemisphere closest to the Earth while the upper trace is from the conjugate 
hemisphere. The simulated echoes in Figure 8B can be compared with the observed traces in 
Figure 2A and show the excellent qualitative agreement between the corresponding local 
hemisphere (A) and conjugate-hemisphere (B) traces. 

4. SIGNAL GUIDANCE MECHANISM 

As mentioned earlier, Reinisch et al., [2001] suggested another possible mechanism to produce 
field-aligned echoes. In that scenario, field-aligned propagating radio signals would naturally 
follow the magnetic field line when the density (refractive index) gradient direction is parallel to 
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the geomagnetic field line. This condition may be likely at high latitudes where the field lines 
and density gradients are more or less in the radial direction as discussed in Fung et al. [2003], 
Inversions of multiple plasmagrams of field-aligned echoes during passages at high latitudes 
have allowed Nsumei et al [2003] to determine the polar cap densities under different 
geomagnetic conditions. As also pointed out in Fung et al [2003], however, the density gradients 
near the equator are mostly oblique or perpendicular to the magnetic field directions. This would 
make the density gradient unlikely to be the responsible guiding mechanism for guiding signals 
to and from the conjugate hemisphere. 

Figure 9 shows quantitatively the angles between the model density gradients and the model 
local geomagnetic field lines as a function of magnetic latitude along the L =2, 2.5, 3, 3.5 and 4 
dipole field lines. As shown in Figure 9, the electron density gradients in the plasmasphere (L < 

4) become highly oblique to the geomagnetic field lines as the magnetic latitudes are less than 
~35°, indicating that signal guidance by field-aligned density-gradient is an unlikely mechanism 
that would lead to the discrete echoes (especially the ones from conjugate hemisphere) often 
observed by RPI in the plasmasphere. On the other hand, the density gradients (thus refractive 
index gradients) at sufficiently high L values, e.g., over the polar cap, are more easily to become 
field-aligned. 

The diffuse echoes, as shown in Figure 2a, that often accompany the field-aligned echoes 
observed in RPI plasmagrams can be viewed as corroborating evidence for the presence of field- 
aligned density irregularities. Carpenter et ah, [2002] and Fung et al [2003] attributed the range 
spreading of the diffuse echoes at least in part to coherent scattering from cross-field electron 
density irregularities of < 10% of the background density that range in size from 200 m to greater 
than a few km. The ray tracing calculations in this paper clearly demonstrates that under these 
conditions we would expect to observe ducted echoes. 
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5. DISCUSSION 

Because of the rather extensive histoiy in investigations of ionospheric ducts, it may be of 
interest to compare the duct density model employed here [equation (6)] and those use in 
ionospheric studies [e.g., James, 2000]. Equation (6) gives the cross-field electron density profile 
of a magnetospheric density duct along a dipole geomagnetic field line. The density structure 
along the field line length is assumed to be in diffusive equilibrium. The analytic prescription is 
used to facilitate the ray-tracing calculations presented earlier. As shown by equation (6), both 
the total density perturbation (SNNde) and the duct width c determine the cross-field density 
gradient of the duct, which controls the amount of refraction imposed by the duct on a given 
propagating wave (f, \|/). Although a constant a is used in the present study, the changing Nde 
along the magnetic field actually provides a relatively realistic treatment of the physical 
situation. An extension of the model is to incorporate a a that varies with distances along the 
field, but we expect those results to be essentially the same as here. 

In a recent study of ionospheric ducts using bi-static radio observation, James [2000] used a 
cylindrical duct model with exponentially decreasing electron density with altitude. The cross- 
sectional duct density profile is modeled by a form factor given by {1- cos[7t(r/rD) m ]}, where r is 
the radial distance from the duct center, rp is the cylindrical duct radius, and m specifies the 
steepness distribution of the cross-field density gradients. Unlike equation (6), this model is valid 
only for modeling field-aligned density ducts in limited altitude range in the low-altitude polar 
ionosphere with fixed scale height and essentially vertical (radial) field-line geometry. 

The Gaussian cross-field density profile in equation (6) is known to be a “well-behave” function 
with no sharp turns that seems to lead readily to ducting of waves with wave-normal directions 
less than -30° from the magnetic field orientation (see Figure 3). On the other hand, the James 
[2000] model tends to require relatively steep duct walls (m = 4 and m = 6) to trap waves even at 
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a few degrees from the magnetic field. The model predicts no effective wave ducting for gentle 
sloping (m = 1) of density or obliquely propagating waves (e.g., \|/ > few degrees). 

We should also examine the conditions valid for ducting by FAI in some details. As shown in 
Figure 10, radio signals is expected to be trapped by field-aligned density structures in which the 
refractive index w(p) encountered by the wave is a decreasing function of the distance p 
(measured normal to B) from the duct center. This occurs for high-frequency R-X and L-O mode 
waves when SAWde < 0, and for whistlers when 8 AWde > 0. As the wave propagates from a 
region of higher n(p,) (e.g., at points A, C, E, and G) to a region of lower «( P 2 ) (at points D, H, 

B, and F), refraction of the wave would turn the wave propagation directions from more oblique 
to more parallel to B. Total internal refraction or trapping of the wave occurs when the wave 
normal angle satisfies \\r < (90° - <j> c ), where <f> c is the critical incidence angle at which the 
refracted wave becomes parallel to the density duct or the background magnetic field. For proper 
ducting, the field-aligned density duct must also satisfy following conditions: 

(1) Reflecting surfaces (shaded regions in Figure 10) along the length of the duct must be 
sufficiently smooth that they are all specular surfaces, i.e., the amplitudes of any surface 
roughness must be small compared to the wavelengths of the ducted waves. Otherwise, 
refraction and scattering by surface roughness could easily prevent ducting over long 
distances. 

(2) The field-aligned density scale length LN/A^ong the duct must be large compared to the 
duct width a, which is the cross-field scale over which N e may vary by a few percent. 
Furthermore, we must also have X < <r(« Lm) for wave propagation. Thus for a given duct 
width a, there is a minimum wave frequency^, ~ c/a below which no waves can be ducted. 

Based on conditions (1) and (2) above, we can estimate the smoothness and width of the 
density duct that may be responsible for the discrete echoes shown in Figure 2a. It is noticed that 
the discrete traces associated with the local hemisphere (trace A) and that reflected in the 
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conjugate hemisphere (trace B) have low-frequency limits at 269.8 kHz and 280.6 kHz, 
respectively. The low-frequency limit of trace A is observed with very small virtual range (~ 0.3 
Re) and therefore must be near the local R-X cutoff frequency at the spacecraft with a free-space 
wavelength X ~ 1.1 km, which is also the limit of surface roughness along the density duct in the 
local hemisphere. The slightly higher low-frequency limit (280.6 kHz) of the conjugate trace (B), 
however, suggests that signal guidance by the conjugate portion of the density duct is not so 
effective at lower frequencies perhaps due to surface roughness of the order of 1 km scale. Both 
traces A and B have a maximum frequency /Max of 369.2 kHz (with X ~ 0.81 km), thus we have 
fudfex. ~ 1.37. This corresponds to a minimum duct size a ~ 107.Max> 81 km (for 1% depletion, 
see Figure 6). 

5. CONCLUSIONS 

The ray tracing calculations presented in this paper show that field-aligned density depletions in 
the plasmasphere could lead to ducting of high-frequency extraordinary mode waves (Figures 3 
and 4). As.shown in Figures 5 and 6, proper guidance of wave at a given frequency 0 m3x//rx) 
and wave normal angle (0 < \|/ < 30°) requires the duct width to increase with increasing density 
depletion (SNNde)- Ducting of nearly field-aligned propagating waves (\j/ ~ 0°) can occur within 
density depletions as small as 1% and less than 10 wavelengths wide (Figure 6). When launched 
at high altitudes, only nearly field-aligned propagating rays can be trapped and guided. 

Conjugate field-aligned echoes observed by RPI are most likely due to waveguide ducting with 
relatively smooth ducts, as prescribed by conditions (1) and (2) above, maintained along much of 
the plasmaspheric field lines. 
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Figure Captions 

Figure 1. A schematic of ducting of a high-frequency (f»fuh) electromagnetic wave (green 
arrows) by field-aligned electron density depletions. 

Figure 2. Panel A shows a typical RPI plasmagram with field-aligned and diffuse echoes 
labeled. Panel B illustrates the directions of the observed echoes relative to the spacecraft 
position depicted by the red dot. As shown in panel A, the diffuse echoes are spread in range and 
are therefore reflecting from a series of irregularities as shown in panel C. 

Figure 3. R-X mode rays at 200-kHz with y/ = 0-360°in 10° intervals, launched at 20° magnetic 
latitude and L = 3.5, 5 N/N D e = -15%, c = 30 km. Only rays with initial yr< 30° or with^> 150° 
are guided by the model density duct. 

Figure 4. Ducted and non-ducted ray at f= 120 kHz, Sy/= 8°, 0° magnetic latitude, L = 

3.5, 3N/Nde= -5%, <r= 150km {60X ) having different wave normal angles (/from 0° to 72° at 
8° intervals launched from the magnetic equator toward the northern hemisphere. Rays with 
large wave normal angles (y/> 30°) are not ducted. 

Figure 5. Maximum ducted frequency 3 of waves (with different initial wave normal 

angles y/) as a function of normalized duct width r| for a 2% N e depletion (panel A) and a 10% 

N e depletion (panel B). Larger density depletions can support sounder waves at progressively 
higher frequencies. 

Figure 6. Variation of maximum ducted frequency P =7max//rx as a function of duct width for 
0° propagation in ducts with different N e depletions. 
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Figure 7. Variation of virtual ranges of echoes at f/ftot =1.1 and different initial wave normal 
angles guided by a field-aligned density duct with rj = 50 and SNNm < 5%. 

Figure 8. Simulated plasmagram traces when sounding from the equator (panel A) and from 
20°magnetic latitude (panel B) in a field-aligned density duct with 5% depletion and r\ = 50. 

Figure 9. A plot of the angle between the model density gradient (VN) and the geomagnetic 
field line (B) along the L =2, 2.5, 3, 3.5, and 4 field lines as a function of magnetic latitude. 

Figure 10. A schematic showing the guidance of signals by total internal reflections along a 
density duct, where the refractive index n as seen by the wave along the duct is a decreasing 
function of the distance from duct center. As the wave propagates along the duct, refraction of 
the wave changes the wave normal angle *F from a finite value near the duct center to 0° (for 
parallel propagation) or 180° (for anti-parallel propagation) near the edge of the duct. 
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